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ABSTRACT 

We report the possible discovery of multiple gravitationally-lensed images of QSOs 
with angular separations on arcminute scales. The QSOs were selected from the com- 
pleted 2dF QSO Survey as having redshifts and optical colours which are statistically 
consistent. In this paper we present higher-quality optical spectra of the candidates 
and discuss the likelihood of their genuinely being systems lensed by massive clusters 
of galaxies. From a comparison of the spectra it appears that up to six pairs of QSOs 
may be lensed multiple images, although the true number may be less than that and 
further observations should be undertaken to amass more evidence and to detect the 
lensing clusters. Two of the candidates may be associated with low redshift clusters 
of galaxies. 
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1 INTRODUCTION 

Strong gravitational lensing of background objects by mas- 
sive clusters of galaxies is expected to produce multiple im- 
ages with wide separations, > 30", and the multiple im- 
ages of background lensed galaxies have long been known 
with radii of giant arc systems up to 70 " (e.g. iKneib et alJ 
ll99dlSahu et al]ll998tlSmail et al.lll995tl . However, despite 
a nu mber of searches, primarily in large-area radio surveys 
(e.g. IMaoz fc Rbdll99a iMaoz et alJll997t lOfek et alJl200lL 
2002: lPhillips et al.ll2001allblh widely-separated multiple im- 
ages of QSOs have not yet been found. This is attributable 
to the low probability of strong lensing along any given line 
of sight coupled with the relatively low numbers of known 
QSOs. The wi dest separation QSO l ensed systems have sepa- 
rations ~ 7" ijKochanek et al.ll2002lh Other wide-separation 
multiple QSO systems found to date are more likely to be 
"binary QSOs" (that is, a pair of active black holes cohab- 
iting the same parent h ost galaxy) rather tha n multiple im- 
ages of the same QSO llKochanek et al.lll999l) . 

We have searched a catalogue of 22,163 QSOs that com- 
prises the recently completed 2dF QSO survey ("the 2QZ 
survey") for candidate wide-separation lensed QSO systems 
and have obtained higher signal-to-noise spectra in an at- 
tempt to eliminate obvious non-lensed systems and to iden- 
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tify systems that are likely to be genuine lensed systems. 
After an incomplete survey of the most likely candidates, 
there remain six systems which, with varying degrees of con- 
fidence, are likely to be lensed. 



2 THE SELECTION OF CANDIDATE LENSED 
SYSTEMS FROM THE 2QZ SURVEY 

The 2QZ survey covers 740 deg 2 in two strips of sky, one 
centred on the South Galactic Pole region and the other 
on the celestial equator in the region of the N orth Galactic 
Pole. The survey photometry is described bv ISmith et alJ 
( 2002) and the initial rel e ase of t he first 10,000 QSOs is de- 
scribed by ICroom et al ] fcOOlal) . In constructing the 2QZ 
survey, candidate QSOs were first selected from it, bj and r 
photographic photometry of morphologically-stellar objects, 
and objects whose colours lay bluewards of the main stel- 
lar locus were selected for low resolution spectroscopy at 
the two-degree field (2dF) facility at the Anglo- Australian 
Telescope. 

The initial spectra obtained were of sufficient quality 
that reliable identifications of QSOs and measurements of 
their redshifts could be made, with a n r.m.s. error in red - 
shift measurement of typically 0.0035 l|Croom et aLlEoOlal h 
although it should be noted that the redshift error for an in- 
dividual object depends on its brightness and emission-line 
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equivalent widths as well as observing conditions. The abso- 
lute accuracy of the photometry in each band is estimated 
as ~ 0.17 magnitudes at the faint survey limit ijSmith et alJ 
2002). The relative photometric accuracy for objects mea- 
sured close together on the same pho tographic plate is better 
than this, however: .Mitchell. |198S) found that, for objects 
more than two magnitudes brighter than the completeness 
limit of the photographic plate, relative photometry accu- 
rate to ~ 0.05 magnitudes was typical, with that value ris- 
ing rapidly to values > 0.1 magnitudes as the plate limit is 
approached. 

The resulting catalogue of photographic photometry 
and spectroscopically confirmed QSOs was then searched 
for pairs or multiple systems whose colours were consistent 
given the large photometric uncertainties and whose red- 
shifts were consistent given the typical redshift uncertainty. 
To determine the colour consistency, it was assumed that 
two QSOs, if lensed, should have the same it, bj and r mag- 
nitudes apart from a magnitude offset arising from differen- 
tial magnification of the images in a lensed system, and thus 
assuming negligible effects due to variations in QSO colours 
with time (see section 0J. The statistic used as the measure 
of colour consistency was 

3 

vIbr = ^(Am, - Ami) 2 , 

i = l 

where Arrii is the difference in apparent magnitude of the 
two QSOs measured in waveband i and Am, is the average 
magnitude difference for that QSO pair. Redshift consis- 
tency was measured by the difference, |Az|2qz in 2QZ cat- 
alogue redshift of two QSOs. Candidate multiple systems 
were initially selected if they satisfied the criteria a UBR < 
0.2, |Az|2Qz < 0.01 and angular separation A9 < 6'. The 
colour criterion is close to the expected scatter in this statis- 
tic, given the photometric errors at the survey limit, and 
hence we expect the initial colour selection to select only 
~ 50% of the lensed systems. The redshift criterion is a fac- 
tor two larger than the expected scatter in Az|2Qz and this 
criteria should allow the majority of lensed systems to be 
selected. The stringent colour criterion was chosen in this 
initial study in order to maximise the probability of find- 
ing lensed systems, and hence further observations should 
be carried out to establish the true completeness of the se- 
lected sample. 

A further check was made on the consistency of their 
radio properties by searching for detections in the NVSS 
JCondon et alJll99Sft . One candidate pair with clearly dis- 
crepant radio detections was removed from the sample of 
candidates. The NVSS properties of two remaining pairs are 
discussed further in section 3. 

Candidates with 18.25 < b.j < 20.85 were selected from 
12307 QSOs in the SGP region and from 9856 QSOs in the 
NGP region of the 2QZ survey, and 38 pairs of QSOs in 
total met these criteria. Of the QSOs with confirmed red- 
shifts, there were no candidate systems selected with more 
than two members, although some of the pairs do have addi- 
tional ufar-selected companions for which 2QZ spectroscopic 
observations either were not carried out or did not lead to 
a positive identification. 

Our expectation is that the majority of the selected 
pairs will not prove to be lensed systems. Test samples 



were also constructed in attempt to assess the possible con- 
tamination. First, QSOs pairs that are more widely sepa- 
rated, but still within one degree of each other and mea- 
sured on the same photographic plates, were selected with 

Az|2Qz < 0.02 and the distribution of (Tubr was com- 
pared with the a UB R distribution of the candidate pairs. No 
statistically-significant difference was found, although this 
is probably a reflection of the large errors in Oubr arising 
from the large photographic photometry errors. A second 
test sample comprising QSO pairs with cj UB r > 0.5 and 

A2I2QZ < 0.01 was selected and the distribution of angular 
separations, normalised by the total number of pairs selected 
with AO < 10', was compared with the candidate lens sam- 
ple selected to have a UBR < 0.2. This time there was an 
excess of pairs on arcminute scales in the sample of candi- 
date lensed systems consistent with a minority of the pairs 
being lensed systems, but because of the small numbers of 
pairs selected the test was not statistically significant. Hence 
at this stage there is no evidence either for or against the hy- 
pothesis that the selected sample contains a subset of lensed 
systems. Further evidence needs to be collected on all mem- 
bers of the sample in order to establish its completeness and 
contamination. 

A key part of that evidence is the follow-up spectroscopy 
described in the next section. Because of limitations of ob- 
serving time and observing conditions, only a subset of the 
candidates have had such observations, and only that subset 
is described in the remainder of this paper. For prioritising 
the observations, the pairs were ranked in \Az\ and also in 
(Tubr- a composite rank was produced by multiplying these 
ranks and the top candidates in each region were observed. 
In the SGP region, 10 pairs were observed, but the spectro- 
scopic follow-up is particularly sparse in the NGP region, 
where only for one pair have spectra been obtained of suf- 
ficient quality to allow some assessment of whether or not 
it may be lensed. At this stage the sample reported here 
cannot be regarded as either being complete or free from 
contamination by non-lensed pairs of QSOs. We postpone 
any discussion of the completeness of the sample to later 
work. We anticipate that in future the candidate QSO pairs 
in the NGP region will also have the benefit of more ac- 
cura te, five-band, colours from the Sloan Digital Sky Sur- 
vey llStoughton et al]l2002l) . which should lead both to sam- 
ples of candidates with lower contamination and with better 
quantified photometric uncertainties. 

In the following sections each pair is designated by the 
IAU naming convention applied to the mean position of the 
pair, in the FK5 (equinox J2000) coordinate system. 



3 FURTHER OPTICAL SPECTROSCOPY 

Spectra were obtained of the candidate pairs at the William 
Herschel 4.2-m. telescope (WHT) on 2002 February 5-6 and 
at the Anglo- Australian 3.9-m. telescope (AAT) on 2002 
September 6-10. At the WHT, the ISIS long-slit spectro- 
graph with 300B grating and EEV detector was used with 
on-chip binning to provide spectra from the blue arm with 
1.7A pixeP 1 dispersion. Observing conditions were poor, 
with seeing of about 2" for the spectra reported here, re- 
sulting in spectra with resolution about 8A. Data were also 
obtained with the ISIS red arm but were of insufficient qual- 
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Table 1. Observed QSO pairs. The columns are: mean right ascension of the pair (equinox 2000), mean declination, mean 2QZ 
survey redshift, colour difference statistic, 2QZ redshift difference and angular separation, followed by the right ascension, declination 
and 2QZ catalogue apparent b magnitude for each member of the pair. 
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ity to be useful. Observations at the AAT were carried out in 
generally good conditions with seeing in the range 0.7—1.6". 
The 300B grating on the RGO long-slit spectrograph with 
EEV detector provided a dispersion of 1.6 A pixel -1 and res- 
olution of typically 6A(depending on the choice of slit width, 
which was varied according to the seeing at the time). 

The typical integration time on each QSO varied be- 
tween one hour and four hours depending on observing con- 
ditions and brightness of the QSO. In most cases the QSOs 
in each pair were sufficiently far apart that they could not be 
observed simultaneously. Instead, integrations on each QSO 
were split into intervals of 30 minutes duration, and no more 
than three integrations were carried out sequentially on one 
QSO without then observing the other member of the pair. 
By this process we could be sure that the spectra were taken 
under sufficiently similar conditions that any differences be- 
tween them are likely to be genuine. For those QSOs ob- 
served in this way the slit position angle was arranged to 
be at the parallactic angle. In two cases, J1435+0008 and 
J0135— 3021, the pair members were close together and were 
observed simultaneously with the slit oriented at the appro- 
priate angle. In the case of J1435+0008 the pair could be 
observed at hour angles such that the slit position angle 
was within 20 degrees of the parallactic angle. In the case 
of J0135— 3021 this was not possible, but the pair could be 
observed within 30 degrees of the zenith where the effects of 
differential refraction should be small. 

Spectra of all the pairs observed are presented in FigQ 
In order to improve the signal-to-noise ratio per pixel the 
spectra have been smoothed with a top-hat function to a 
resolution specified in the figure caption for each pair. Each 
panel shows spectra of the two members of each pair super- 
imposed. The spectra have been scaled by the total value 
of the observed (uncalibrated) count rate so that they may 
be compared. Each panel also shows the difference between 
the scaled spectra, for comparison. The spectra have been 
bias-frame subtracted and flatfielded in order to remove in- 
strumental features from the spectra, but have not had any 
flux calibration applied, in order to facilitate evaluation of 
any differences. The noise in the spectra is dominated by 
photon noise from the sky background, and the noise in each 
wavelength pixel is more uniform with wavelength in uncal- 
ibrated spectra. Absorption features arising from the terres- 
trial atmosphere remain in the spectra, most prominently 



at A7600A, and should be ignored when comparing spectra. 
The ratios of the total count rate between the spectra are 
given in Tabled 

One expectation of the "lensed QSO" hypothesis is 
that the two images should have the same redshift. As the 
pairs were selected to have 2QZ survey redshift differences 
|Az|2qz < 0.01 we would, in the absence of QSO clustering, 
expect the distribution of WHT or AAT redshift differences 
to be uniformly distributed within that range. QSO cluster- 
ing has a scale lengt h in redshift of about 0.003 at z ~ 1.5 
JCroom et al]l2001rl . and if the QSOs in each pair are in- 
dividual objects we would expect them also to have redshift 
differences arising from their peculiar velocities, so the "sep- 
arate QSOs" hypothesis predicts that even in the presence of 
QSO clustering there should still be a broad distribution of 
redshift differences when measured from the WHT or AAT 
spectra. 

To test this, the redshift differences between each pair 
of WHT or AAT spectra have been measured by cross- 
correlation, about subtracting the continuum from the un- 
calibrated spectra and avoiding the atmospheric absorption 
feature at A7600A. Uncertainties in the redshift differences 
have been estimated assuming the noise in the spectra is 
dominated by uncorrelated photon shot noise in each pixel. 
A Monte-Carlo method was adopted in which Gaussian noise 
of the same amplitude as the noise in each spectrum was 
added to smoothed versions of each spectrum to create a 
number of synthetic spectra with the same spectral fea- 
tures and noise properties as the actual data. The synthetic 
spectra were then cross-correlated with each other and the 
quoted uncertainty in the redshift difference is the rms of the 
values derived from the synthetic spectra. These uncertain- 
ties are lower limits to the true uncertainty as they contain 
no allowance for systematic differences and noise that may 
be correlated between pixels. However, the spread in values 
of the measured redshift differences is consistent with these 
uncertainty estimates, indicating that photon shot noise is 
the dominant source of uncertainty. Remarkably, we find six 
of the pairs to have velocity differences less than 200 km s" 1 
(Tableland Fig.[5J. This is a strong indication that a num- 
ber of these pairs are lensed systems. All separations for 
these pairs are sufficiently large that, if they were two dif- 
ferent QSOs, the linear separations would be greater than 
1 Mpc and they must occupy separate host galaxies. QSO 
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Figure 1. Scaled and superimposed spectra of the candidate lensed systems. The differences between the scaled spectra are 
also shown. To help distinguish genuine features from noise, all spectra are shown without any relative flux calibration applied. 
Atmospheric absorption features are still present in the uncalibrated spectra. The flux scale is the number of detector counts 
observed per pixel in individual exposures (typically of duration 1800 s) of one of the QSOs, with the count values of the second 
QSO scaled to match, upper left WHT ISIS blue-arm spectra of J1435+0008 smoothed to a resolution of 12A. upper right 
AAT RGO spectrograph spectra of J0150— 2921 smoothed to a resolution of 12A. lower left AAT RGO spectrograph spectra of 
J0135— 3021 smoothed to a resolution of 15A. lower right AAT RGO spectrograph spectra of J0215— 2906 smoothed to a resolution 
of IK A 
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Figure 1. (continued): left AAT RGO spectrograph spectra of J0223-2949 smoothed to a resolution of 15A. right AAT RGO 
spectrograph spectra of J2228— 2857 smoothed to a resolution of 15A. 



clustering would cause a peak in the redshift distribution at 
Az| = 0, but broadened not only by the measurement error 
but also by the relative peculiar velocities of the QSO host 
galaxies. Peculiar velocities at high redshift are expected to 
be lower than in low-redshift galaxy surveys, but simulations 
still predict pairw ise velocity dispersions a ~ 350kms _1 
jZhao et alJ 12003) . significantly larger than would be in- 
ferred from Fig. [S] 

Although the redshift difference is a quantitative mea- 
sure of similarity between the spectra, we should also look 
for other differences in the spectra. We now give a brief 
evaluation of the comparison of the spectra of each pair. 
The discussion is ordered by an approximate degree of con- 
fidence in the likelihood that the pair is a lensed QSO. This 
evaluation is based not only on redshift difference but also 
on similarity of continuum shapes, emission line equivalent 
widths, line ratios, velocity widths and profiles. This com- 
parison is necessarily somewhat subjective. It is never possi- 
ble to prove that a pair of spectra are multiple images of the 
same lensed QSO, and in this case we have to rely on build- 
ing up a body of evidence which can then be interpreted 
as either supporting the lensed hypothesis or providing ev- 
idence against it. The spectra we have obtained only cover 
the range of accessible optical wavelengths and vary widely 
in their signal-to-noise. The spectra of a QSO pair viewed 
with high signal-to-noise may appear to be discrepant, such 
as J2250— 3105, whereas such a discrepancy may not have 
been apparent in lower signal-to-noise spectra. The only reli- 
able solution to this effect is to obtain higher quality spectra 
and other evidence for those pairs with poor signal-to-noise. 



With these factors in mind, any evaluation at this stage 
cannot be regarded as conclusive, but nonetheless should 
help us to identify those QSOs that are most likely to be 
lensed, that can then be the subject of further observations 
to amass more evidence. The candidates we have observed 
were selected a priori to have similar optical colours (albeit 
with a large uncertainty) and redshifts, so we should not be 
too surprised to find that almost all of the observed pairs 
do indeed have very similar overall continuum shapes in the 
observations presented here. Instead, we must appeal to the 
wide range of detail in QSO spectra which we would not ex- 
pect to match when comparing the spectra of independent 
QSOs. 

3.1 J1435+0008 

This pair was the only one presented here that was observed 
at the WHT. Other pairs were attempted but only inad- 
equate spectra were obtained. Observing conditions were 
poor, with poor atmospheric transparency and poor seeing. 
The ISIS red arm spectra contain too little signal, and we 
show here a comparison of the blue-arm spectra only. The 
dichroic between the red and blue arms of ISIS causes a 
sharp fall in response at A > 5000A, with some oscillations 
in that response in this wavelength range. This is a good 
lensed QSO candidate, with a small redshift difference and 
an unusual ratio of Ly-a and NV which is common to both 
spectra. 

This QSO pair also appears i n the Sloan Digital Sk y 
Survey (SDSS) Early Data Release dStouehton et al]l200a\ 
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Figure 1. (continued): upper left AAT RGO spectrograph spectra of J2235— 2717 smoothed to a resolution of 15A. upper right 
AAT RGO spectrograph spectra of J2210— 3204 smoothed to a resolution of 12A. lower left AAT RGO spectrograph spectra of 
J2203— 3133 smoothed to a resolution of 12A. lower right AAT RGO spectrograph spectra of J2250— 3105 smoothed to a resolution 
of 12A. 
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Figure 1. (continued): Uncalibrated AAT RGO spectrograph 
spectra of J2329— 2957 smoothed to a resolution of 12A. 



Table 2. Observed redshift differences and flux ratios between 
the WHT/AAT spectra. The flux ratio is defined as the ratio of 
the flux from the second listed QSO in Table 1 to the flux from 
the first listed QSO of each pair. 
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Figure 2. The histogram of WHT/AAT redshift differences, 
\Az\, listed in Tabled The original selection limit based on 2QZ 
redshifts is shown, as is the redshift difference expected for a pair 
at the mean sample redshift of z = 1.8 with a relative peculiar 
velocity difference of 200 km s — 1 



lensed-QSO hypothesis: if these are images of the same QSO 
then, at redshifts near the QSO, the Ly-a forest absorp- 
tion lines should be identical in both spectra. If the pair 
were not a single lensed QSO, then the QSO separation 
would be ~ 1 Mpc and the Ly-a lines would not be ex- 
pected to be common between the QSOs: sight lines that are 
separated on such scales when comparing unrelated QSOs 
show significan t differe nces in their Ly-a forest absorption 
jD'Odorico et alll998Tl . Forest lines at redshifts significantly 
different from the QSO may nonetheless differ between the 
images in either hypothesis because the light rays for each 
image would have diverged to a significant extent. In the 
case of J 1435+0008 the ultra-violet part of the spectrum is 
noisy, but nonetheless it does appear to have possible ab- 
sorption features in common between the spectra. 
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Table E| gives the SDSS photometry for this pair, together 
with the differences in apparent magnitudes and the associ- 
ated error. It can be seen that the colours of the two QSOs 
are consistent with being the same to substantially higher 
accuracy than was possible with the 2QZ photographic pho- 
tometry, although there does seem to be some some evidence 
for a systematic difference in colour with waveband. This 
may be taken as evidence against the system being lensed, 
or it may be accounted for either by spectral variability or 
by differential extinction between the light paths. 

We also note the possibility of an elegant test of the 



3.2 J0150 2921 

This is possibly one of the most enigmatic pairs of spectra 
obtained. The spectra have good signal-to-noise ratio and 
appear identical in every respect apart from one gross fea- 
ture: the presence of a broad absorption line system in one 
of the spectra. For reasons described in section [I] we do not 
believe that this should be taken as evidence against the 
system being a lensed QSO, and because of the otherwise 
striking similarity between the spectra in these observations 
with good signal-to-noise we regard this as being one of the 
systems most likely to be a single lensed QSO. 

This system also has a number of bright galaxies in 
the vicinity, of which three have h ad redshifts measur ed by 
the 2dF Galaxy Redshift Survey (IColless et al.ll2003l) . The 
positions and redshifts of the two QSOs and of the three 
galaxies are shown in Fig. [3] This system is discussed further 
in section 5. 

We also note here that although neither of the QSOs is 
detected in the NVSS survey there are four detected sources 
of radio emission extending over some arcminutes approxi- 
mately symmetrically placed with respect to the centre of 
the system, and approximately perpendicular to it. It may 
be that the radio emission is the substructure of a single ex- 
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tended double radio source associated with one of the galax- 
ies in the field. 

3.3 J0135 3021 

There are some small differences visible in the ultraviolet 
parts of the spectra of this pair which lead to the somewhat 
high value of redshift difference given in Table H However, 
these spectra were taken simultaneously, with the slit ro- 
tated to the position angle defined by the coordinates of 
the two QSOs, and the observations presented here were 
made with the pair within about thirty degrees of the zenith. 
Slight errors in aligning the slit with the pair could result 
in a small differential loss of light owing to atmospheric re- 
fraction which may contribute to the small differences at 
A < 4500A. At this stage we still regard this pair as possi- 
bly being a lensed QSO, although further observations are 
needed. 

3.4 J0215 2906 

This pair has very similar spectra with a redshift difference 
of zero to high accuracy. There are some residual features 
present in the difference of the scaled spectra which may ar- 
gue against the lensed QSO hypothesis but which may sim- 
ply be observational noise or may reflect low-level emission 
variations (see also section QJ. 

3.5 J0223-2949 

This pair has spectra identical within the noise except in the 
CIV line, where there is a residual difference whose maxi- 
mum is about one-third the maximum of the emission-line. 
The other emission lines and the detailed continuum shape 
are extremely similar, however, and the redshift difference 
is consistent with zero at the 2-sigma level. We regard this 
as being a possible lensed system that requires further ob- 
servations. 

This system also has a number of bright galaxies in the 
vicinity (Fig. of which ten have ha d redshifts measure d 
by the 2dF Galaxy Redshift Survey dColless et al]l2003t) . 
This system is discussed further in section 5. 

There is a 8.6 mjy NVSS (ICondon et al]|l998l) radio 
source 54" to the north of the northern QSO which is not 
reproduced near the southern QSO. However, that radio 
source is coincident with the bright peculiar-shaped galaxy 
at z = 0.1147 and it is unlikely that the radio source is 
related to the QSO. 

There are also weaker NVSS detections at the NVSS 
detection limit, one a 2.8 mjy source 45" from the south- 
ern QSO, the other a 3.5 mjy source 24" from the northern 
QSO. Both the reality of these weak radio detections and the 
possibility of their being associated with the QSOs needs in- 
vestigation with a higher resolution, higher sensitivity radio 
map. 

3.6 J2228-2857 

This pair has very similar spectra but does display some 
significant differences in the Ly-a/NV and CIV lines. This 
could be a lensed system if we are prepared to accept 



the possibility of time variations in emission-line strength 
and/or QSO intrinsic absorption features (see section^J. 



3.7 J2235 2717 

Although the redshift difference for this pair is consistent 
with being zero, the QSOs do have significantly different 
emission-line ratios, and we believe this pair is rather un- 
likely to be a lensed system. 



3.8 J2210-3204 

This pair exhibits significant differences in line strength be- 
tween the two spectra and is also unlikely to be a lensed 
QSO. 



3.9 J2203 3133 

This pair has quite different continuum shapes and a large 
redshift difference and is unlikely to be a lensed QSO. 



3.10 J2250 3105 

This pair has very similar spectra, but there does appear to 
be a continuum difference between the ultraviolet parts of 
the spectra and a significant redshift difference. We consider 
this unlikely to be a lensed pair. 



3.11 J2329 2957 

This pair exhibits significant differences in emission line 
equivalent widths and ratios and a significant redshift dif- 
ference, making it also unlikely to be a lensed QSO. 



4 TIME DELAY AND THE FORMATION OF 
BROAD ABSORPTION FEATURES 

At least two of the pairs have spectra which match very 
well their redshifts, continuum slopes, emission-line equiv- 
alent widths and profiles but which do nonetheless show 
some significant differences between the two QSOs. One, 
J0150— 2921, has spectra which appear identical apart from 
the existence of CIV broad absorption in one QSO and not 
the other. Such an observation has previously been used 
to argue agai nst QSO pairs be ing lensed versions of the 
same object iGreen et alJ I2002T) but the recent discovery 
that b road abso rption lines can form on timescales of ten 
years (Ma 2002) changes that conclusion. 

In fact, large separation lenses could have a significant 
time delay between the lensed images. The actual time delay 
is strongly dependent on the lens mass distribution and on 
the symmetry of the light paths through the lens potential. 
As we have so far detected a possible lens in no more than 
one of these systems, and as the mass distributions in any 
lenses are unknown, we cannot yet calculate a predicted time 
delay. But in principle we should expect a delay which could 
be as large as the time delay calculated from the maximum 
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possible path length difference in a simple geometry (and 
ignoring the time dilation term) 

< Yc (D - Dlens) 

where 6 is the angular separation of the images, D is the 
observer-QSO angular diameter distance and Dlens is the 
observer-lens angular diameter distance, assuming a flat uni- 
verse. For an angular separation of 200" this corresponds to 
approximately 5000 years. The closeness of the match of 
the spectra of J0150— 2921 indicates that in fact the broad 
absorption features have indeed formed in one QSO within 
this timescale. The available time is certainly much greate r 
than the 10 years formation time observed bv iMal |2002). 
One explanation for such rapid formation of absorption fea- 
tures is the ejection of remnant material from solar-mass 
objects being destroyed in the vicinity of the black hole. As 
the crossing time of the broad-line region of such material 
is of order 10 years this seems a feasible explanation. An al- 
ternative explanation for both t he differenc e in our spectra 
and also the variability seen by IMal i2002T) might be that 
the size of a broad-absorption line gas cloud is only of order 
10 -3 of the distance of the cloud from the central source of 
ultraviolet radiation. 

If this conclusion is correct, it also reopens the possi- 
bility that the 7.3" pair, Q2345+007, is a lensed QSO and 
not a binary QSO. lGreen et all 1120021) have noted that the 
X-ray properties differ, but also note that one QSO shows 
a broad absorption line system. If the difference in X-ray 
properties is caused by X-ray absorption associated with 
the optical broad absorption line feature, as suggested by 
iGreen et alJ i2002f) , then formation of an absorption system 
on a timescale of order 5 years could explain the observed 
results. 

Two other QSO pairs exhibit some differences in their 
broad emission lines, and again it is possible th at these also 
arise from long-timescale spectral variability. Sma ll et alJ 
il997l) have studied spectral variability specifically with this 
question in mind, and find that emission-line variations of 
this magnitude are common on timescales of a few years. 

Time variation may have the unwanted effect of causing 
the sample to be incomplete, if colours or spectral features 
vary significantly so that either candidates are not initially 
selected, or if they are subsequently rejected in a spectral 
comparison. Given the potentially long time-delays between 
light paths in systems of such wide separation, it is difficult 
to allow for this effect as no measurements of variability in 
unlensed systems can be made on those timescales. In the 
long term the best approach would be to select candidates 
on the basis of observables that should be time-independent, 
such as redshift difference between multiple images, and to 
apply independent tests of the lensing hypothesis such as the 
Ly-a forest comparison discussed above or looking for char- 
acteristic features of lensing such as giant arcs associated 
with a QSO's host galaxy. In this paper we do not address 
this question further, preferring instead to view this paper 
as providing evidence for the existence of extreme systems 
which may then be used as evidence to justify more complete 
surveys. 



Table 3. SDSS Early Data Release photometry for the QSO 
pair J1435+0008, showing the apparent magnitudes in the SDSS 
wavebands, their associated photometric errors, and the magni- 
tude differences between the QSOs and the errors in those quan- 
tities. 



QSO 


u 


9 


r 


i 


z 


J143508+000844 
± 


21.37 
0.09 


20.72 
0.03 


20.57 
0.03 


20.64 
0.05 


20.37 
0.13 


J143506+000901 
± 


20.49 
0.05 


19.92 
0.02 


19.87 
0.02 


19.98 
0.03 


19.74 
0.08 


Am 
± 


0.88 
0.10 


0.80 
0.04 


0.70 
0.04 


0.66 
0.06 


0.63 
0.15 



5 IDENTIFICATION OF THE LENSES 

If these systems of paired QSOs are indeed lensed QSOs, 
then the lenses should be identifiable as clusters of galax- 
ies. As described in section 2, one of the possible lensed 
systems, J0223 — 2906, appears associated with a galaxy 
cluster at z = 0.11. Ten galaxies in the vicinity have red- 
shifts, measured by the 2dF Galaxy Redshift Survey, that 
are close to this value, with a velocity dispersion of about 
600 km sec -1 (Fig.|3J. Such a velocity dispersion would how- 
ever produce an image splitting of only about 20" for an 
isothermal sphere lens, so if this cluster is the lens then ei- 
ther an anisotropic or a more centrally concentrated mass 
distribution is needed. The total mass required to produce 
the observed image splitting is ~ 6 x 10 14 MQ. There is also 
a weak ROSAT-detected X-ray source close to the position 
of the two bright galaxies which lie along the line joining the 
two QSOs, although the origin of the emission is not known. 
At present it is not clear if sufficient mass is present in this 
cluster to produce the observed image splitting. One other 
system, J0150— 2921, may also have a low-redshift cluster 
associated with it, again based on redshift measurements 
from the 2dF Galaxy Redshift Survey, although in this case 
only three redshifts are currently known (Fig. [3J . 

In order to estimate the expected distribution 
of lens re dshifts we co ns ider the calculations pre- 
sented by lLopes fc Miller! i2003T) in which the lens- 
ing clusters are assumed to have spherically-symmetric 
iNavarro. Frenk fc White! lll997tl mass profiles, and t o have 
an abundance predicted by the ISheth fc Tormenl (Il999t) 
evolving mass fu nction. Here we assum e a cold dark matter 
power spectrum teardeen et alJll98rj) modified to include 
an approximate correction for baryons JSugivamalli995Tl m 
a ACDM cosmology with Qm = 0.3 and Q,a = 0.7 and a 
power-spectrum normalisation given by as — 0.9. In this 
model, it is found that the distribution of lens redshifts is 
not sensitive to the choice of cosmological parameters, al- 
though the abundance of lensing systems is. Ma g nifica tion 
bias is calculated as discussed bv lLopes fc Milled fcOOSt . 

Within this model, it is clear that a lensing cluster is 
expect to have 2lens < 1, and for the widest separations the 
model predicts that any lens should have zlens < 0.5. How- 
ever, the total probability predicted for 200" lensed systems 
is so low that such systems are expected only if there are 
significant departures from the assumptions of the model, 
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Figure 3. ESO sky survey images of 10' regions around J0150— 2921 and J0223— 2949 showing the locations of both QSOs and 
the redshifts of the galaxies in the vicinity measured by the 2dF Galaxy Rcdshift Survey. The orientation and image scale are also 
shown. 
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Figure 4. The a posteriori probability distribution of redshift, 
•ZLENSi f° r a lensing cluster producing multiple images of a QSO 
at z = 2 with separations either > 200" or > 30", assuming the 
ACDM cosmology and cluster lens model described in the text. 



in which case the lens redshift distribution would need to be 
recalculated. 



6 COMPARISON WITH MODELS 

At this stage it cannot be regarded as conclusive that these 
candidate systems are indeed lensed, and we also do not 
yet have an accurate assessment of the completeness of the 
lensing search. Nonetheless it is interesting to compare the 
res ults presented here wit h the expectations of l ensing mod- 
els. iLi fc OstrikeJ <2002|) and lLopes fc Milieil (|2003h have 
calculated the expected distribution of angular separations 
arising from strong lensing by massive clusters of galax- 
ies in the standard cosmological framework that assumes 
structure grows hierarchically from initial Gaussian pertur- 
bations in a universe dominated by cold dark matter and 
dark energy. With cosmological parameters as determined 



from the WMAP observations JSpergel et alJl2005t the ex- 
pected number of systems with A6 > 10" in the 2QZ survey 
of 2 x 10 4 QSOs is about 0.5 fcopes fc Milled I2003D . so it 
is possible, although not highly likely, that the 2QZ sur- 
vey contains wide-separation lensed systems. However, the 
distribution is expected to be dominated by angular separa- 
tions towards the lower end of the range, and in this stan- 
dard model th e expected number of systems with A9 > 30" 
is only ~ 0.1 l|Lopes fc Millerl|2003F ). One possible explana- 
tion is that the widest-separation systems presented here 
are not in fact lensed, despite the similarity of the redshifts 
and spectra. A second is that the values of cosmological 
and cluster model parameters are not as assumed in calcu- 
lating the above numbers: key parameters are as, w, Qm 
and also parameters of the model such as the dependence 
of cluster halo concentration on mass, redshift and cosmol- 
ogy ancitb_e_£o^niic scatter in halo concentration (for details 
see lLopes fc Miller! I2003T) . Or it may be that the models 
do not yet comprise the full complexity required to model 
accurately the statistics: possible deficiencies in clude the as- 
sumptio n of the | N avarro 1 Frenk fc White! l)l997!) mass profile 
(see also lQguri et al.ll2003l) and the assumption of a smooth 
distribution of mass. Other departures from the standard 
cosmological model, such as the existence of non-Gaussian 
fluctuations, could also modify the lensing probabilities. 

For these reasons we believe it is imperative to test carry 
out observational tests of the candidate lensed systems pre- 
sented here in order to test these models, rather than using 
the predictions of the models to argue against the interpre- 
tation of these systems as being lensed. 



7 COMPARISON WITH PREVIOUS 
OBSERVATIONS 

Previous searches have been undertaken for lensed QSOs 
with comparable separations, and it is interesting to inves- 
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tigate whether the detection of any systems in 2QZ would 
be i nconsistent with t h ose se arches. 

lOfek et al.l i200ll l2002f) searched a sample of size ~ 
8000 radio-loud QSOs without finding any evidence for mul- 
tiple images i n the separation range 5" < AO < 30". 
iPhillips et aflfeOOlhft examined a complete sample of 11,670 
flat-spectrum radio sources for possible multiple images sep- 
arated in angle by 6" < A9 < 15" and found no lensed 
systems. An incomplete sample of ~ 15, 000 radio sources 
yielded one possible lensed system. The first step in com- 
paring these results with those from the 2QZ survey is to 
compare sample selection. 

The range of angular separations probed by the radio 
surveys is generally smaller than considered here, although 
had any 2QZ pairs with A9 > 10" been lens candidates 
they should have been identified as such. For comparison 
purposes we shall use the cluster-lensing model to compare 
surveys covering different ranges in separation, although we 
have argued above that the overall lensing probabilities pre- 
dicted by such models may be too low to explain the exis- 
tence of more than one lensed sy stem in the 2QZ survey. The 
model of lLopes fc M iller (2003) with parameters Om = 0.3, 
SIa = 0.7, (78 = 0.9 and zqso = 1.5 predicts relative lens- 
ing probabilities of 1, 1.7 and 1.3 for the separation ranges 
6 — 15", 5 — 30" and > 10" respectively, assuming the same 
magnification bias in each case. If the models are indeed 
deficient in some way, or if the values of parameters are dif- 
ferent from that assumed, these relative probabilities will be 
different, but they suffice for the discussion in this section. 

In fact, because of the double radio a nd optical selection 
of the FIRST survey dOfek et al.ll200lT) the magnification 
bias could be a factor 5 smaller for that survey compared 
with the JVAS/CLASS or 2QZ surveys (lOfek et al.ll2002t 
IPhillips et all 1200 lri lLopes fc Milled 12003ft . so theoverall 
expected relative lensing probabilities in the three surveys 
are 0.2, 1.7 and 1.3 for FIRST, JVAS/CLASS and 2QZ re- 
spectively. 

We also need to make some assumption about sample 
completeness and to multiply the probabilities by the sam- 
ple sizes. Here we shall adopt a conservative view of the 
JVAS/CLASS result and say that out of 15,000 QSOs no 
lensed systems were found, and we adopt an assumed com- 
pleteness of 50% for the 2QZ survey. At this stage this figure 
is illustrative only. With these assumptions, the predicted 
relative numbers of systems in each survey are 0.1, 1.6 and 
1.0 for the FIRST, JVAS/CLASS and 2QZ surveys respec- 
tively. 

Now suppose the entire dataset contains one lensed sys- 
tem: the a posteriori probability of that system being found 
in the 2QZ survey is 37%. However, if there are more sys- 
tems in the 2QZ survey, the probability of there being none 
in the other surveys falls to 13, 5, 2 % for 2, 3 or 4 2QZ 
systems. Thus, if the surveys are indeed comparable in this 
way, the existence of two lensed systems within the 2QZ sur- 
vey would not be inconsistent with the other surveys, but a 
larger number could be. 

Finally, however, we note that another possibly sig- 
nificant factor in calculating the expected lensing prob- 
abilitie s is the redshift d i stribu ti on of the source p opu- 
lation. IWambseanss et"all J2003ft . lOeuri et all J2003ft and 
lLopes i^Mill^Jll2^0^H iave demonstrated that cluster lens- 
ing probability strongly increases with the mean redshift of 



the source population, and also depends to some extent on 
the distribution of source redshifts about the mean. The 
redshift distribution is not well-known for the CLASS and 
FIRST surveys, but appear s likely to have a median value 
in the range 1.0 < z < 1.4 llOfek et al.ll2002l : IPhillips et alJ 
l2001allbh . The 2QZ median redshift is 1.5. Lensing proba- 
bility increases by a factor about 3 as the source r edshift 
is increased from z — 1 to z = 1.5 jLopes fc Milled l2003ll . 
so if the typical redshift of sources in the radio surveys is 
even slightly lower than that of QSOs in the 2QZ survey 
this would also lower the relative probability of discovering 
lensed systems in the radio surveys. 

Our overall conclusion from this comparison is that, ob- 
servationally, the results from the 2QZ and radio surveys 
may not be incompatible, but clearly further work needs to 
be done both to establish more accurately the redshift dis- 
tribution of the JVAS/CLASS sources and to establish the 
reality of the 2QZ candidate systems. 



8 CONCLUSIONS 

Evidence has been presented for the discovery of the first 
lensed QSOs with multiple images separated on arcminute 
scales. If these are genuine lensed systems then we should ex- 
pect the associated lenses to be observable as massive clus- 
ters of galaxies. Further evidence needs to be collected to 
establish that these are indeed lensed, however. One such 
piece of evidence that is independent of assumptions of the 
nature of the lens would be the existence of common Ly- 
a forest absorption in candidate lensed QSO spectra. The 
existence of giant arcs associated with a system would also 
provide strong evidence for cluster lensing, especially if the 
host galaxies of the QSOs could be detected and shown to be 
distorted by lensing. Detection of the lens itself, either by di- 
rect imaging or weak-lensing, would also provide convincing 
evidence. 

Consideration of the expected distribution of angular 
separations indicates that, although a number of the iden- 
tified systems may be lensed, a number of those with the 
greatest separations may not be. Confirmation or discount- 
ing of the widest separated systems is particularly important 
as any such systems would challenge our current understand- 
ing of lensing by galaxy clusters. 

Among the examples found, one candidate lensed sys- 
tem has been found to display broad CIV absorption in one 
image and not the other, which could be due to the for- 
mation of such absorption systems on timescales of tens of 
years. Further searches for wide-separation lenses may be 
partially incomplete because such absorption systems would 
affect optical colours, spectra and the flux in other wave- 
bands, including possibly the X-ray band. Broad-line region 
variability may also render search methods based on colour 
and/or spectroscopic comparison partially incomplete un- 
less generous selection criteria are allowed. Finally, it may 
be that optical searches are biased against lenses with sepa- 
rations less than a few tens of arcsec because of obscuration 
in galaxies associated with the cluster lens. 
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